AFRL-RY-WP-TP-20 13-0 158 


METAMATERIALS,  SPATIAL  DISPERSION,  AND 
NEGATIVE  REFRACTION  (PREPRINT) 


John  Derov 

Electro-Optic  Components  Branch 
Aerospace  Components  &  Subsystems  Division 

Teresa  O’Donnell 

Air  Force  Life  Cycle  Materials  Command 
Hanscom,  AFB 

Beverly  Turchinetz 
Analog  Devices 


AUGUST  2013 

Interim 


Approved  for  public  release;  distribution  unlimited. 

See  additional  restrictions  described  on  inside  pages 

STINFO  COPY 


AIR  FORCE  RESEARCH  LABORATORY 
SENSORS  DIRECTORATE 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7320 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 

_ Technical  Paper  Preprint 


1.  REPORT  DATE  (DD-MM-YY) 

_ August  20 1 3 _ 


4.  TITLE  AND  SUBTITLE 

METAMATERIALS,  SPATIAL  DISPERSION,  AND  NEGATIVE  REFRACTION 
(PREPRINT) 


6.  AUTHOR(S) 

John  Derov  (AFRL/RYDP) 

Teresa  O’Donnell  (Hanscom,  AFB) 
Beverly  Turchinetz  (Analog  Devices) 


3.  DATES  COVERED  (From  -  To) 

1  October  2011  -  1  June  2013 


5a.  CONTRACT  NUMBER 

In-house _ 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

61102F _ 

5d.  PROJECT  NUMBER 

3001 _ 

5e.  TASK  NUMBER 

12 


5f.  WORK  UNIT  NUMBER 

Y044 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Electro-Optic  Components  Branch  Hanscom  AFB  REPORT  NUMBER 

Aerospace  Components  &  Subsystems  Division  Analog  Devices  AFRL-RY-WP-TP-20 13-0 158 

Air  Force  Research  Laboratory,  Sensors  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7320 

Air  Force  Materiel  Command,  United  States  Air  Force 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSORING/MONITORING 

Air  Force  Research  Laboratory  AGENCY  ACRONYM(S) 

Sensors  Directorate  AFRL/RYDP _ 

Wright-Patterson  Air  Force  Base,  OH  45433-7320  H-  sponsoring/monitoring 

Air  Force  Materiel  Command  AGENCY  REP0RT  NUMBER<S> 

United  States  Air  Force  AFRL-RY-WP-TP-20 13 -01 58 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


13.  SUPPLEMENTARY  NOTES 

PAO  Case  Number  88ABW-2013-3465.  Clearance  Date  6  August  2013.  Report  contains  color. 


14.  ABSTRACT 

Spatial  dispersion,  anisotropy,  and  negative  refraction  with  negative  group  velocity  have  been  observed  in  metamaterials. 

Spatial  dispersion  occurs  when  permittivity  e(ca,  k)  and  permeability  p(to,  k)  are  dependent  on  frequency,  o,  and  wave  vector,  k, 
related  to  a  fixed  periodic  potential  along  a  lattice.  The  potential  along  the  lattice  results  from  an  exciton,  resonance,  or  atomic 
polarization  forming  an  electric  dipole  moment.  The  dipole  moments  produce  a  spatially  dispersive  longitudinal  and  transverse 
wave  allowing  an  incident,  longitudinal  and  transverse  to  exist  in  the  metamaterial.  The  resonance  of  Split-Ring-Resonators 
(SRRs)  in  SRR-Wire-Post  (WP)  medium  provides  the  dipole  moment  needed  for  spatial  dispersion  and  p  <  0  for  negative 
refraction.  The  WP  is  designed  to  have  e  <  0  at  the  resonant  frequency  of  the  SRR  providing  simultaneous  e  <  0  and  p  <  0  for 
negative  refraction.  We  ran  parametric  simulations  changing  the  geometric  parameters  of  the  SRR  and  WP  and  used  the 
simulated  results  to  demonstrate  the  coexistence  of  spatial  dispersion  and  negative  refraction  in  a  SRR-WP  medium.  The  results 
show  the  interaction  of  spatial  dispersion  and  negative  refraction  in  SRR-WP  media.  We  show  this  by  controlling  the 
geometrical  parameters  of  the  SRR  and  WP  conditions  can  be  met  to  differentiate  the  behavior  between  spatial  dispersion  and 
negative  refraction  in  the  SRR-WP  media. _ 

15.  SUBJECT  TERMS 

metamaterials,  spatial  dispersion,  negative  refraction,  split-ring-resonators 


19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

17.  LIMITATION 
OF  ABSTRACT: 


18.  NUMBER 
OF  PAGES 


John  Derov 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

N/A 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


Metamaterials,  Spatial  Dispersion  and  Negative  Refraction 


John  S.  Derov* 

Air  Force  Research  Laboratory,  Sensors  Directorate,  Wright  Patterson  AFB,  OH  45433 

Teresa  H.  O’Donnell* 

Air  Force  Life  Cycle  Materiels  Command,  Hanscom  AFB,  MA  01731 

Beverly  W.  Turchinetz* 

Analog  Devices,  804  Woburn  St.,  Wilmington,  MA  01887 
(Dated:  July  25,  2013) 

Abstract 

Spatial  dispersion,  anisotropy,  and  negative  refraction  with  negative  group  velocity  have  been 
observed  in  metamaterials.  Spatial  dispersion  occurs  when  permittivity  e(u i,k)  and  permeability 
g(u ;,  k )  are  dependent  on  frequency,  ui,  and  wave  vector,  k,  related  to  a  fixed  periodic  potential  along 
a  lattice.  The  potential  along  the  lattice  results  from  an  exciton,  resonance,  or  atomic  polarization 
forming  an  electric  dipole  moment.  The  dipole  moments  produce  a  spatially  dispersive  longitudinal 
and  transverse  wave  allowing  an  incident,  longitudinal  and  transverse  to  exist  in  the  metamaterial. 
The  resonance  of  Split-Ring-Resonators  (SRRs)  in  SRR-Wire-Post  (WP)  medium  provides  the 
dipole  moment  needed  for  spatial  dispersion  and  g  <  0  for  negative  refraction.  The  WP  is  designed 
to  have  e  <  0  at  the  resonant  frequency  of  the  SRR  providing  simultaneous  e  <  0  and  g  <  0  for 
negative  refraction.  We  ran  parametric  simulations  changing  the  geometric  parameters  of  the  SRR 
and  WP  and  used  the  simulated  results  to  demonstrate  the  coexistence  of  spatial  dispersion  and 
negative  refraction  in  a  SRR-WP  medium.  The  results  show  the  interaction  of  spatial  dispersion 
and  negative  refraction  in  SRR-WP  media.  We  show  this  by  controlling  the  geometrical  parameters 
of  the  SRR  and  WP  conditions  can  be  met  to  differentiate  the  behavior  between  spatial  dispersion 
and  negative  refraction  in  the  SRR-WP  media. 

PACS  numbers:  78.61.Pt 
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I.  INTRODUCTION 


In  Veselago’s  seminal  paper1  on  Left-Handed  material,  he  described  a  material  with  simul¬ 
taneously  negative  permeability  (/i)  and  permittivity  (e)  in  a  left-handed  coordinate  system. 
The  Left-Handed  material  described  by  Veselago  with  a  negative  refractive  index  has  two 
required  conditions,  e  and  ji  must  be  simultaneously  negative  and  it  must  have  a  negative 
group  velocity.  Veselago  also  indicated,  aside  from  a  negative  refractive  material,  materials 
exhibiting  anisotropy  and  spatial  dispersion  have  a  negative  group  velocity  as  well.  In  2000, 
the  first  left-handed  material,  or  metamaterial,  with  simultaneous  negative  e  and  /i ,  was 
demonstrated  using  Split- Ring- Resonators  (SRRs)  and  Wire-Post(WP)2-4.  At  resonance, 
the  SRRs  provides  negative  q5  and  the  WP  provides  negative  e6  producing  negative  refrac¬ 
tion  in  a  metamaterial,  which  has  been  demonstrated  in  the  literature3.  Veselago’s  model 
only  accounted  for  the  frequency  dependence  of  e(ui)  and  in  a  metamaterial.  To  de¬ 
scribe  the  behavior  observed  for  e  and  /i,  a  bi-anisotropic  model1'8  and  spatial  dispersion9  12 
models  were  used.  Spatial  dispersion  has  provided  the  best  explanation  for  the  behavior 
of  e  and  /i  in  a  SRR-WP  metamaterial.  Alu’s  model11,12  used  both  e  and  n  to  determine 
frequency  and  spatially  dependent  behavior  of  a  metamaterial.  Since  spatial  dispersion  and 
negative  refraction  are  dependent  on  the  resonance  of  the  SRRs  and  negative  group  veloc¬ 
ity  in  the  metamaterial,  SRR-WP  media  were  chosen  to  study  the  interaction  between  the 
spatial  dispersion  and  negative  refractive  index  in  metamaterials.  Periodic  SRR-WP  media 
were  used  to  determine  how,  or  if,  spatial  dispersion  can  be  differentiated  from  negative  re¬ 
fraction  in  a  metamaterial.  The  dispersion  in  a  typical  material  e{oj)  and  n(cu)  is  dependent 
on  frequency  uj. 

A  spatially  dispersive  material  occurs  when  e(co,  k)  and  k)  are  dependent  on  fre¬ 
quency  u  and  a  wave  vector  k,  where  k  is  related  to  a  fixed  periodic  potential  along  a  lattice. 
The  periodic  potential  in  a  medium  can  be  produced  by  an  exciton,  resonance,  or  a  polar¬ 
ized  lattice  at  the  molecular  level.  The  exciton,  resonance,  or  polarized  lattice  form  electric 
dipole  moments  at  the  lattice  sites  of  a  medium.  In  our  metamaterial,  the  SRR  resonance 
provides  the  fixed  potential  with  period  a  of  the  SRR-WP  unit  cell,  which  determines  k, 
where  k  =  2n/a  meeting  the  condition  for  spatial  dispersion  in  the  medium.  Since  the 
conditions  for  spatial  dispersion  and  negative  refraction  are  related  to  e  and  /q  we  used  e 
and  q  to  study  the  spatial  dispersion  and  negative  refraction  in  metamaterials. 
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A  diagram  of  the  balanced  SRRs  with  a  WP  used  for  this  study  is  shown  in  Figure  1. 
As  seen  in  Figure  la,  the  balanced  SRRs  consist  of  2  identical  rings  parallel  to  each  other 
with  their  gaps  opposing  each  other.  Opposing  currents  are  produced  in  the  rings,  creating 
symmetric  field  distributions  around  the  SRRs.  The  WP  is  placed  between  the  2  SRRs  for 
symmetry  as  shown.  The  orientation  of  the  linearly  polarized  plane  wave  used  to  excite  the 
SRR-WP  metamaterial  is  also  shown  in  Figure  la.  The  electric  field  must  be  parallel  to  the 
WP  in  order  to  produce  a  negative  e.  Figure  lb  shows  the  unit  cell  with  cubic  symmetry, 
which  was  used  for  all  the  simulated  structures.  The  unit  consists  a  dielectric  0.50  mm 
thick  with  a  relative  e  =  2.2  to  support  the  SRRs  and  WP  where  the  SRRs  reside  on  the 
outer  surface  and  the  WP  is  in  the  middle  of  the  dielectric.  The  rest  of  the  unit  cell  volume 
consists  of  air. 


Plane  Wave 


(a) 


a=b=c=  4.1  mm 


(b) 


FIG.  1.  (a)  SRRs  with  the  WP  between  them  and  the  linear  polarized  plane  wave  used  to  excite 
the  SRRs  resonance,  (b)  SRR-WP  unit  cell  used  for  the  simulations. 


In  earlier  work13,  we  fabricated  SRR-WP  media  and  measured  their  transmittance.  Snell’s 
law  was  used  to  determine  the  refractive  index  for  the  SRR-WP  medium.  The  transmittance 
and  refractive  index  for  the  measured  and  simulated  SRR-WP  medium  and  unit  cell  respec¬ 
tively  is  shown  in  Figure  2.  The  unit  cell  for  the  medium  in  Figure  2  was  simulated  using 
Ansyst’s  High  Frequency  Structure  Simulator  (HFSS)  and  periodic  boundary  conditions. 
Though  the  simulated  results  in  Figure  2  are  for  the  unit  cell  rather  than  the  medium,  com¬ 
paring  the  trend  in  the  frequency  response  between  the  simulated  and  measured  results  show 
enough  agreement  to  use  SRR-WP  unit  cell  simulations  to  demonstrate  spatially  dispersive 
and  negative  refractive  behavior  in  the  SRR-WP  media. 
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FIG.  2.  Measured  and  Simulated  transmission  response  and  real  part  of  the  refractive  index. 

To  look  at  the  spatially  dispersive  behavior  and  negative  refractive  index  of  the  SRR-WP 
structure,  parametric  simulations  were  run  with  varying  geometric  dimensions  of  the  SRR 
and  WP,  holding  the  dimension  of  the  unit  cell  constant.  The  dimensional  changes  of  the 
SRRs  and  WP  used  for  the  simulation  are  in  Table  I. 


TABLE  I.  Dimensions  for  SRRs,  WPs,  and  Unit  Cell 


Dimension 

OD  mm 

ID  mm 

Gap  mm 

Unit  Cell 

SRR  Size 

2.79 

2.54 

0.51 

4.06  mm 

Width  mm 

Length 

WP  Size 

0.75 

ps  oo 

0.25  ,  0.23,  ...0.025 

PS  oo 

SRR  Size 

3.81 

3.56 

0.51 

4.06 

WP  Size 

0.75 

~  oo 

II.  SPATIAL  DISPERSION 


Scattering  parameters  for  the  SRR-WP  unit  cells  were  obtained  from  the  HFSS  simu¬ 
lations.  Using  a  Nicolson-Ross  extraction14,  e  and  p  were  determined  from  the  scattering 
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parameters  of  the  unit  cell.  The  result  of  the  extraction  for  a  unit  cell,  with  a  2.79  mm 
Outer  Diameter  (OD)  SRR  and  0.75  mm  wide  WP  are  shown  in  Figure  3.  The  frequency 
and  spatial  dependence  of  e(cu,  k)  and  h(lu,  k)  are  observed  for  the  unit  cell  in  Figure  3. 


Unit  Cell  with  2.79  mm  OD  SRRs 


Frequency  (GHz) 


FIG.  3.  Spatial  and  frequency  dispersion  e(u;,k),  /i(cj,  k),  e(u;),  and  /x(u ;)  for  the  unit  cell  with  the 
2.79  mm  OD  SRRs  with  0.75  mm  WP. 


To  differentiate  the  spatial  and  frequency  dependence  in  e(cj,  k)  and  /i(cj,  k)  the  frequency 
dependence  of  e(cu)15  and  r(cu)5  were  calculated  for  the  unit  cell.  Equation  1  was  used  to 
calculate  e(co)  for  the  WP,  where  cop  is  the  plasma  frequency  for  the  WP  and  fi(co)  was 
calculated  using  equations  2,  where  uim0  is  the  resonant  frequency  of  SRRs.  To  fit  /i(cu)  to 
/i(cu,  k),  F  and  7  were  varied  until  fi(cu)  matched  um0  and  ump  where  comp  is  the  magnetic 
plasma  frequency. 


rM 


e(w)  =  1  -  uj^/uj2 

(1) 

Fu2/(u2  -  u^o  +  tw) 

(2) 

The  spatial  dispersion  can  be  observed  by  comparing  the  difference  in  the  line  shape  between 
e(cu,  k)  and  e(u>)  and  /d(u>,  k)  and  h(uj)  in  Figure  3.  A  resonance  is  observed  in  the  line  shape 
for  e(co,  k)  that  does  not  occur  in  e(co)  and  the  resonance  in  e(co,  k)  coincides  with  the 
resonance  in  /r(u;,  k),  indicating  there  is  an  interaction  between  the  SRRs  and  WP.  The 
effect  of  the  spatial  dispersion  on  /i  can  be  seen  in  the  deviation  of  the  line  shape  for  /x(c 0, 
k)  from  the  Lorentzian  line  shape  of  //(cu).  To  determine  why  the  resonance  is  occurring  in 
e(uj,  k)  and  how  it  relates  to  the  coupling  between  the  WP  and  SRRs,  a  SRR-WP  unit  cell 
was  simulated  without  the  WP  to  look  at  the  resonant  behavior  of  the  SRRs.  The  results  of 
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2.79  mm  OD  SRRs  without  WP 


FIG.  4.  e(w,  k),  n(ui,  k),  and  /i( cu)  for  the  2.75  OD  SRRs  without  a  WP. 


the  simulation  for  the  SRRs  are  shown  in  Figure  4.  Figure  4  shows  there  is  a  resonance  in  e 
and  e  >  0,  which  coincides  with  the  resonance  in  /i  for  the  SRRs.  Indicating  the  resonance 
observed  in  e  for  the  SRR-WP  in  Figure  3  is  associated  with  the  SRRs.  Since  the  resonance 
for  e  >  0  in  Figure  4  and  e  <  0  for  the  resonance  observed  for  the  SRR-WP  unit  cell  in 
Figure  3  requires  coupling  between  the  SRRs  and  the  WP  for  the  resonance  to  occur  when 
e(cu,  k)<  0.  Figure  4  also  shows  the  spatial  dispersion  is  associated  with  the  SRRs  since 
k)  deviates  from  the  Lorentzian  line  shape  of  /i(ca)  for  the  SRRs  between  com0  and  ump, 
which  correlates  with  the  frequency  range  of  the  resonance  in  e.  The  resonance  in  e(cu,  k) 
of  the  SRR-WP  unit  cell  in  Figure  3  shows  the  WP  is  coupling  with  the  SRRs,  the  effect  of 
the  coupling  between  the  SRRs  and  WP  was  studied  by  changing  the  width  of  the  WP. 

The  width  of  the  WP  was  varied  from  0.75  mm  to  0.025  mm  for  the  2.79  mm  OD  SRR- 
WP.  Figure  5a  and  b  shows  the  real  part  of  e(ca,  k),  k),  e(ca),  and  /i(u;)  for  a  0.25  and 

0.13  mm  WP  width  respectively.  In  Figure  5a  the  line  shape  of  k)  is  in  good  agreement 
with  the  Lorentzian  line  shape  of  /x(c a).  Comparing  the  resonance  of  e(c a,  k)  in  Figure  5a 
and  3,  the  resonance  has  decreased  and  e(cu,  k)  is  approaching  e(u>)  below  uip. 

The  reduction  of  the  intensity  of  the  resonance  for  e(c a,  k)  as  the  WP  width  was  decreased 
and  the  agreement  of  /i(cu,  k)  with  indicates  the  spatial  dispersion  has  been  reduced 

to  a  level  where  it  can  be  considered  negligible.  In  Figure  5b,  the  line  shape  for  /i(u;,  k)  does 
deviate  from  the  Lorentzian  line  shape  of  /j.(co)  and  the  resonance  in  e(u ;,  k)  shows  better 
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8  and  ja  for  SRR-WP  with  0.25  mm  Wide  Post 


3, 

3 

CD 

oi 

oQ 


O' 

0 

01 


Frequency(GHz) 

(a) 


0 

oi 

oQ 


0 

01 


(b) 


FIG.  5.  Real  Part  of  e(ui,  k),  n(ui,  k),  e(w),  and  n{u)  for  a  Unit  Cell  with  2.79  mm  OD  SRRs  with 
a)  0.25  mm  WP,  b)  0.13  mm  WP. 


agreement  with  e(ui)  then  e(oj,  k)  in  Figure  5a.  Though  the  spatial  dispersion  is  negligible 
in  Figure  5b,  value  of  e(uj,  k)  is  approaching  0  between  cumo  and  uimp  as  the  width  of  the  WP 
decreased.  As  e(u,  k)  — »  0  the  value  of  the  negative  refraction  index  -n  — >■  0. 

As  stated,  early  spatial  dispersion  occurs  due  to  the  resonance  in  e(u,  k),  which  creates 
an  electric  dipole  moment  along  the  WP  in  the  unit  cell  producing  the  spatial  dispersion  in 
the  SRR-WP  medium.  In  the  metamaterial  there  is  a  longitudinal  and  a  transverse  wave 
associated  with  the  spatial  dispersion  and  an  incident  wave  propagating  through  the  medium. 
The  incident  wave  propagating  through  the  SRR-WP  medium  excites  the  resonance  in  the 
SRRs,  producing  the  electric  dipole  moments  along  the  WP  and  the  longitudinal  wave  in 
the  medium.  The  period  of  the  longitudinal  wave  in  the  medium  is  given  by  a  the  period 
of  the  SRR-WP  unit  cell  and  provides  the  spatial  dependence  for  e(o ;,  k)  where  k  =  27 r/a. 
The  frequency  of  the  longitudinal  wave  is  the  same  as  the  resonant  frequency  of  the  SRRs 
even  though  the  period  of  the  longitudinal  wave  is  a.  The  incident  and  longitudinal  wave 
propagate  through  the  medium  at  the  same  frequency,  but  at  different  periods.  The  period 

7 


Approved  for  public  release;  distribution  unlimited. 


used  for  the  incident  wave  is  its  free  space  wavelength9  A  =cq/ f  where  cq  is  the  speed  of  light 
in  vacuum  and  /  is  the  frequency  of  the  incident  wave.  As  the  incident  wave  propagates 
through  the  medium,  the  phase  of  the  incident  wave  and  the  phase  of  the  longitudinal  wave 
interact.  This  is  analogous  to  the  way  a  plane  wave  propagates  through  a  crystalline  lattice. 
The  core  wave  functions  of  atoms  along  the  lattice  support  the  phase  of  the  wave  as  it 
propagates  through  the  crystal16  as  shown  in  Figure  6. 


Plane  Wave 


Propagating  Plane  Wave  and  Atomic  Core  Wave  Function 


FIG.  6.  Plane  wave  propagation  along  a  crystalline  lattice  as  illustrated  by  Ziman. 


Consider  a  SRR-WP  unit  cell  in  a  medium  as  the  atoms  in  a  crystal  lattice  and  the 
resonance  of  the  SRR-WP  are  its  core  wave  functions.  The  interaction  between  the  phase  of 
the  incident  wave  and  the  longitudinal  wave  is  shown  in  Figure  7.  The  period  and  resonant 
frequency  SRR-WP  unit  cell  in  Figure  3  was  used  to  determine  the  phase  response  of  the 
medium  in  Figure  7,  making  the  frequency  of  the  incident  wave  11.4  GHz  (A=26.3mm)  and 
the  period  of  the  longitudinal  wave  4.1  mm.  In  Figure  7,  the  maxima  and  minima  for  the 
phase  of  the  incident  wave  do  not  lineup  with  the  maxima  and  minima  of  the  longitudinal 
wave.  The  mismatch  in  the  phase  of  the  incident  wave  with  the  longitudinal  wave  shows 
the  incident  wave  is  not  totally  supported  by  the  longitudinal  wave,  as  the  incident  wave 
propagates  through  the  medium,  which  produces  a  phase  mismatch  between  the  incident 
and  longitudinal  wave. 

To  show  the  propagation  of  an  incident  wave  that  is  supported  by  the  longitudinal  wave 
in  the  medium,  the  4.1  mm  period  of  the  unit  cell  and  the  condition  a/ A  <  1  for  spa¬ 
tial  dispersion9  was  used  to  determine  the  resonant  frequency  of  the  SRRs.  The  resonant 
frequency  determined  for  the  SRRs  was  7.2  GHz(A=41.7nnn).  In  Figure  8,  the  phase  propa¬ 
gation  of  the  7.2  GHz  incident  wave  and  4.1  mm  longitudinal  wave  for  a  SRR-WP  medium. 
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Sinusoidal  Phase  Propation  in  a 
2.79  mm  OD  SRR  and  0.75  mm  WP  Medium 


FIG.  7.  Phase  propagation  of  the  11.4  GHz  incident  wave  and  the  4.1  mm  longitudinal  wave  in  a 
SRR-WP  medium. 

Figure  8,  the  maxima  and  minima  of  the  phase  for  the  incident  wave  does  lineup  to  the 
maxima  and  minima  of  the  phase  of  the  longitudinal.  The  phase  match  between  the  2  waves 
indicates  that  the  phase  of  the  incident  wave  is  supported  by  the  longitudinal  wave,  as  it 
propagates  through  the  medium. 


Sinusoidal  Phase  Propation  in  a 
3.81  mm  OD  SRR  and  0.75  mm  WP  Medium 


FIG.  8.  Phase  propagation  of  the  7.2  GHz  incident  wave  and  the  4.1  mm  longitudinal  wave  in  a 
SRR-WP  medium. 

To  match  the  resonant  frequency  SRRs  to  4.1  mm  period  of  the  unit  cell,  the  outer 
diameter  of  the  SRRs  was  increased  to  3.81  mm.  A  unit  cell  for  the  3.81  mm  OD  SRRs  with 
a  0.75  mm  wide  WP  was  simulated  and  e((u,  k),  k),  e(u>),  and  /i{uj)  were  determined 

from  the  results  of  the  simulation.  The  real  part  of  e(cu,  k),  /i(cu,  k),  e(u>),  and  /i(cu)  for  the 
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unit  cell  is  shown  in  Figure  9.  The  resonance  of  e(uj,  k)  is  significantly  stronger  for  the  7.2 
GHz  unit  cell  then  e(u>,  k)  for  the  11.4  GHz  unit  cell  in  Figure  9. 


8  &  \A  for  3.81  mm  OD  SRRs 
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FIG.  9.  The  real  part  of  e(w,  k)  and  k)  with  the  e(u,  k)  resonance  strongly  couple  to  the  WP. 


III.  NEGATIVE  REFRACTIVE  INDEX 

We  have  discussed  spatial  dispersion  in  SRR-WP  media  and  found  that  the  spatial  disper¬ 
sion  can  be  controlled  using  the  geometric  dimensions  of  the  SRRs  and  WP.  We  also  found 
that  negative  refraction  and  spatial  dispersion  are  not  independent  of  each  other  since  the 
negative  refractive  index  and  spatial  dispersion  occurs  within  the  frequencies  of  resonance 
for  f i .  The  negative  refractive  index  occurs  for  frequencies  where  e(u ;,  k)  and  k)  are 
simultaneously  negative,  the  refractive  index  n—y/ql,  n  was  determined  using  e(co,  k)  and 
h(oj,  k)when  e(ui,  k)  and  k)  simultaneously  negative  (-)n  was  assigned  to  the  refractive 
index.  Figure  10  shows  the  refractive  index  for  the  2.79  mm  OD  unit  cells  where  the  WP 
width  was  changed  between  0.75  and  0.05  mm.  In  Figure  3,  the  frequencies,  which  refractive 
index  is  negative,  are  between  ujm0  and  comp  except  for  the  0.05  mm  WP. 

Taking  a  closer  look  at  how  the  coupling  between  the  SRRs  and  WP  is  effected  by  the 
change  in  the  WP  width,  we  will  look  at  the  frequencies  where  e  and  /i  are  simultaneously 
negative  with  respect  to  the  WP  width.  Using  2.79  mm  OD  SRR  unit  cells  and  their  values 
of  u)mo,  Ump ,  and  ujp,  the  values  of  um0  were  plotted  as  a  function  of  the  WP  width.  Figure 
11  shows  UJmo,  tOmp ,  and  ojp  as  a  function  of  the  WP  width.  The  frequencies  when  fi  and  e 
are  simultaneously  negative  are  shown  in  Figure  11.  In  Figure  11,  p  <  0  between  (um0  and 
ujmp  and  for  c op  >  comp  e  <  0  over  all  the  frequenneies  for  p  <  0.  If  ujp  <  ojmp  and  up  >  (jJmo 
then  e  <  0  over  only  part  of  the  frequency  range  between  umo  and  c omp  as  seen  in  Figure 
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FIG.  10.  Refractive  Index  for  the  2.71  mm  OD  SRRs  with  WP  0.75,  0.25,  0.13,  and  0.05  mm 


10  for  the  0.05  mm  wide  WP.  Figure  11  also  shows  the  frequency  varies  for  ujm0  and  ump  as 
the  WP  width  changes,  but  Ac op  =  c omp  -  um0  almost  constant  and  Ac op  is  approximately 
equal  to  u>mp  -  cum0  of  the  SRRs  without  the  WPs.  Showing  the  change  in  the  width  of  WP 
causes  um0  and  ump  for  the  SRRs  to  shift  in  frequency  relative  to  (jm0  and  c omp  for  the  SRRs 
without  the  WP.  The  interaction  between  the  SRRs  and  the  WP  changes  when  WP  =  0.25 
mm.  For  WP  widths  >  0.25  mm  cumo  and  c omp  are  greater  than  their  values  for  the  SRRs 
without  a  WP  and  for  WP  widths  <  0.25mm  cumo  and  c omp  are  less  than  their  values  for  the 
SRRs  without  the  WP. 

Plasma  &  Resonant  CO  2.79  mm  OD  SRRs 


FIG.  11.  SRR  ump  and  ujm o  for  the  2.79  mm  OD  SRRs  and  WP  cjp  for  widths  from  0.75  to  0.025 
mm  and  the  displaced  ujp  for  the  0.25  mm  WP  width. 


To  look  at  the  effect  of  the  SRRs  on  the  WP,  ojp  was  calculated  for  the  WP  widths 
without  the  SRRs6.  Equation  3  was  used  to  determine  the  c op  for  the  WP  in  the  unit  cell, 
where  neff  is  the  effective  charge  density  for  WP,  e  is  the  electron  charge,  er  is  the  effective 
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relative  permittivity,  for  the  unit  cell  used  in  the  simulations  to  account  for  the  dielectric 
used  to  support  SRRs  and  WP  along  with  the  with  the  volume  of  air  in  the  unit  cell,  eo  is 
the  permittivity  of  free  space,  and  meff  is  the  effective  electron  mass  for  the  wire. 

„2 


u2p  = 


n. 


ffe 


6r€()lTlef f 


(3) 


The  neff  is  determined  by  the  spatial  fraction  that  the  WP  occupies  in  the  unit  cell  and  is 
given  in  Equation  4  where  n  is  the  electron  charge  density  in  the  WP,  L  is  the  width  of  the 
wire,  W  is  the  thickness  of  the  wire,  and  a  is  the  period  the  unit  cell. 


LW 

Ueff  ~  n~^2~ 


(4) 


The  effective  meff  is  given  by  equation  5, 

p0  LWe2n 

meff  = - - ln(a/r)  (5) 

where  po  is  the  permeability  of  free  space  and  r  =  \J~LW~Jt:  .  Substituting  equations  4  and 
5  into  3,  then  up  is  given  by  equation  6  for  WP  in  the  unit  cell. 


Co 


2vr 


(jJp 


a  y  er  In  ( a/r ) 


(6) 


where  c0  is  the  speed  of  light  in  vacuum.  Figure  12  shows  the  calculated  up  for  the  WP 
without  the  SRR  and  u>p  for  the  simulated  SRR-WP  as  a  function  of  WP  width.  Looking  at 
the  curves  in  Figure  12  for  ujp,  all  the  plasma  frequencies  calculated  using  the  WP  widths 
nicely  fall  on  the  curve  as  would  be  expected. 


Post  Width  (mm) 


FIG.  12.  Plasma  frequency  vs.  post  width  for  SRR-WP  unit  cell  for  post  widths  from  0.75  mm  to 
0.025  mm  and  calculated  for  post  wide  L=0.75  to  0.025  mm,  W=  0.05mm,  er=l.ll. 
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The  curve  of  the  SRRs-WP  unit  cells  shows  a  response  similar  to  the  calculated  plasma 
frequencies  for  the  WP  widths  form  0.025  to  0.23  mm  indicating  the  WPs  are  behaving  more 
like  their  calculated  counterparts,  until  the  width  of  the  WP  is  0.25  mm.  For  the  WP  widths 
0.25  and  0.75  mm  up  shows  a  stronger  interaction  with  the  SRRs.  The  up  of  0.25  and  0.75 
mm  WP  widths  for  the  SRR-WP  unit  cells  show  a  larger  deviation  from  the  calculated  up 
for  the  WP  widths.  Both  the  0.25  and  0.75  mm  WP  widths  show  a  strong  interaction  with 
the  SRRs,  bnt  the  0.25  mm  WP  width  shows  the  strongest  interaction  with  the  SRRs.  In 
Figure  11  the  0.25  mm  up  for  the  0.25  mm  WP  width  approaches  ump  for  the  SRRs  where 
the  value  of  up  for  the  0.75  mm  WP  width  is  closer  to  its  calculated  u>p. 

In  Figure  4  we  showed  that  the  resonance  in  e  was  associated  with  resonance  of  p  for 
the  SRRs  and  the  SRRs  were  responsible  for  the  spatial  dispersion  in  the  SRR-WP  meta- 
material.  In  Figure  5a  we  showed  that  the  interaction  between  the  0.25  mm  WP  virtually 
removed  the  spatial  dependence  in  p  by  its  agreement  with  p(o;).  The  disassociation  up  for 
the  0.25  mm  wide  WP,  from  its  calculated  behavior  in  Figure  12  and  the  result  from  Figure 
5a  indicates  the  behavior  of  e  and  p  is  localized  in  the  unit  cell.  Since  p  does  not  have  a 
spatial  dependence,  it  allows  e  and  p  to  be  treated  independently  in  Maxwell’s  equations, 
indicating  a  medium  made  from  the,  2.79  mm  OD  and  0.25  mm  WP,  units  cells  can  be 
homogeneous.  Treating  p  as  just  a  frequency  dependent  parameter  provides  a  way  to  dis¬ 
tinguish  the  negative  refraction  in  a  metamaterial  even  though  e  even  though  e  is  spatially 
dependent. 


IV.  SUMMARY  AND  CONCLUSION 

As  discussed,  the  waves  produced  by  spatial  dispersion  in  the  metamaterial  have  a  longi¬ 
tudinal  and  transverse  mode.  A  plane  wave  incident  on  a  SRR-WP  metamaterial  excites  the 
resonance  in  the  SRRs  and  couples  to  the  WP  creating  an  electric  dipole  moment,  producing 
a  longitudinal  wave  in  the  metamaterial.  The  coupling  between  SRRs  and  WP  allows  energy 
from  the  incident  wave  to  be  transferred  to  the  longitudinal  wave  and  the  2  waves  interact  as 
they  propagate  through  the  metamaterial.  When  the  wavelength  of  the  resonant  frequency 
of  the  SRRs  and  the  period  of  the  unit  cell  coincide,  the  phase  of  the  longitudinal  and  the 
incident  wave  will  match,  increasing  the  coupling  between  the  SRR  and  WP  producing  a 
stronger  electric  dipole  moment  along  the  WP.  The  increased  strength  of  the  electric  dipole 
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moment  in  the  unit  cell  shows  the  excitation  of  the  SRRs  by  incident  wave  is  transferring 
its  frequency  domain  energy  to  the  longitudinal  spatial  domain  wave  via  /j  by  the  resonance 
in  e. 

We  also  showed  the  association  of  negative  refraction  with  spatial  dispersion  in  a  meta¬ 
material  since  the  both  the  negative  refractive  index  and  the  spatial  dispersion  depend  on 
the  resonance  of  the  unit  cell.  Varying  the  width  of  the  WP  allowed  us  to  control  the  spatial 
dispersion  in  the  SRR-WP  medium  and  by  matching  the  resonant  frequency  of  the  SRR- 
WP  with  the  period  of  the  unit  cell  the  spatial  dispersion  increased  in  a  SRR-WP  medium. 
Though  the  spatial  dispersion  can  be  maximized  and  minimized  in  the  SRR-WP  medium, 
the  negative  refraction  and  spatial  dispersion  could  not  explicitly  be  separated.  By  varying 
the  width  of  WP  in  the  unit  cell  we  determined  that  a  0.25  mm  wide  WP  minimized  the  spa¬ 
tial  behavior  of  the  SRR-WP  unit  cell  with  the  2.79  OD  mm  SRRs.  The  coupling  between 
SRRs  and  the  0.25  mm  WP  in  the  unit  cell,  virtually  eliminates  the  spatial  dependence  in  /i 
making  the  negative  refractive  index  distinguishable  from  the  spatial  dispersion.  Concluding 
the  spatial  dispersion  and  negative  refraction  in  a  SRR-WP  cannot  be  explicitly  separated, 
but  negative  refraction  can  be  distinguished  from  the  spatial  dispersion. 
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